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We unveil the general features of the phase diagram for any gauge theory with fermions transform- 
ing according to distinct representations of the underlying gauge group, at the four-loop order. We 
classify and analyze the zeros of the perturbative beta function and discover the existence of a rich 
phase diagram. The anomalous dimension of the fermion masses, at the infrared stable fixed point, 
are presented. We show that the infrared fixed point, and associated anomalous dimension, are well 
described by the all-orders beta function for any theory. We also argue the possible existence, to all 
orders, of a nontrivial ultraviolet fixed point for gauge theories at large number of flavors. 
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Determining the phase structure of generic gauge the- 
ories of fundamental interactions is crucial in order to be 
able to select relevant extensions of the standard model 
of particle interactions [ 1 J . In particular we are interested 
in elucidating the physics of non-Abelian gauge theories 
as function of the number of flavors, colors and matter 
representation. 

To gain a quantitative analytic understanding of the 
phase structure of different gauge theories we investi- 
gate the zeros of the perturbative beta function to the 
maximum known order and for one of the zeros also the 
limit of large number of flavors to all-orders. 



I. ZEROLOGY 

We consider the perturbative expression of the beta 
function and the fermion mass anomalous dimension 
for a generic gauge theory with only fermionic matter in 
the MS scheme to four loops which was derived in f2H3l: 



da 



dlnu 2 

X 

S3 



= j8(a) = -pV - fro 3 - pV - pV + 0(a 6 ) , (1) 



dh\m 
din y 2 



y(a) 



= y a + y x a 2 + y 2 a 3 + y 3 a 4 + 0(a 5 ) , (2) 



where m = m(fi 2 ) is the renormalized (running) fermion 
mass and /.< is the renormalization point in the MS scheme 
and a = a/ An = g^/lSn 2 where g = g(ii 2 ) is the renor- 
malized coupling constant of the theory. 

The explicit expression of the coefficients above are 
reported in the appendix [A| f or completeness. Note also 
that the beta function is gauge independent, order by 
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order in perturbation theory |2 j. The same also holds for 
the anomalous dimension of the fermion mass y. 

Here we investigate the structure of the zeros of the 
four-loops beta function for any matter representation 
and gauge group. Interestingly we find a universal clas- 
sification of the behavior of the zeros as function of the 
number of flavors rif. 

The first observation is that due to the fact that the 
beta function is a polynomial of degree five in a there are 
five complex zeros. Since one can factor out a 2 , the beta 
function will always have a double zero at the origin. The 
other three zeros determine the interesting properties of 
the theory, to this loop order. In the following we will 
focus on these three zeros which can be either all real or 
one real and two complex. 

To elucidate the landscape of possible topologies we 
plot the real nontrivial zeros as function of the number 
of flavors normalized to the one above which asymptotic 
freedom is lost (rif) in Fig.JT] Solid black lines correspond 
to the location of the ultraviolet (UV) zeros while red- 
ones to the infrared (IR) stable fixed points. The shaded 
areas denote the regions where the |3 function is positive. 
We have rescaled the vertical axis using the function 
a* = 2 arctan (5a) / n, mapping [-co, +oo] into the interval 
[-1,1]. 

The best way to read these figures is to imagine a 
straight vertical line corresponding to a fixed value of 
rif. The intersection of this line with the solid curve de- 
termines the number of the zeros, the color of the curves 
the type of zeros (if red is IR and if black is UV), and fi- 
nally the corresponding horizontal value is the coupling 
location. We term the landscape of the zeros the zerology. 

We investigate also the negative values of a since this is 
the most natural mathematical setting. In fact, the prop- 
erties of the pure Yang-Mills theory at negative a have 
already been studied on the lattice by Li and Meurice in 
|4 |. There the authors have shown interesting relations 
between the positive and negative regions of a. 

It is possible to identify four distinct topologies able to 
fully represent the entire zerology for any gauge grou 
and matter representation which are reported in 
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(a) First kind of topology. 
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(b) Second kind of topology. 




(c) Third kind of topology. 




(d) Fourth kind of topology. 



FIG. 1: The four different topologies displayed above classify the entire zerology landscape. We show, in each plot, 
the regions of positive (gray) and negative (white) values of the beta function for different gauge theories. The 
solid lines, per each figure, are the locations of the zeros of the beta functions. The lines of UV fixed points are in 

2 

black while the IR ones in red. We have defined a* - — arctan (5a). The vertical dashed red-lines correspond to the 

71 

location where one zero approaches infinity. 



We start by summarizing a number of general features 
that we have identified: 

• At small number of flavors there is only a negative 
ultraviolet zero. 

• At around and above nj we observe the existence 
of three zeros, two ultraviolets and one infrared. 
The infrared one, near Hj, is the Banks-Zaks J5| 
point. Above rip the IR fixed point is now at a 
negative values of a and at a new critical number 
of flavors collides with the UV fixed point zero at 
a negative value of the coupling, forming a double 
zero. At this point the beta function is positive for 
any negative alpha. 

• At very large number of flavors the U V fixed point, 
for positive values of alpha, always exists and ap- 
proaches zero asymptotically as n~ 2 ^ 3 . The explicit 



derivation is provided in Sec. Ill 



• By increasing n f from zero there is always a critical 
number of flavors above which an IR fixed point 
emerges for positive a. 

The distinguishing feature of different topologies is how 
the zeros merge or disappear as function of nj. 

The topology A (Fig. la I is characterized by the fact 
that the zeros always remain at finite values of the cou- 
pling. This means that when a zero disappear it has to 
annihilate with another one. This happens at two dis- 
tinct locations. One at a positive value of the coupling 
and the other at a negative one occurring after asymp- 
totic freedom is lost. 

In the topology B, represented in Fig. lb as for the 
previous case, we still observe the merging of the IR 
and UV zeros at two different number of flavors. In this 
case, however, there is a region in the number of flavors, 
where the UV fixed point located at positive couplings 
reaches infinity at finite n t and appears on the negative 
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Rep. | Top. A | Top. B | Top. C | Top. D 
SWN) 



FUND 


N = 2,3 


4<N< 11 


N> 12 




ADJ 








N> 2 


2-SYM 








N> 2 


2-ASY 


N = 3,4,5 


N = 6,7 


8 < N < 26 


N>27 


SO(N) 


FUND 




N< 6 


N = 5 


N = 3,4 


ADJ 








N> 3 


2-SYM 








N> 3 


SP(2N) 


FUND 


N= 1,2 


3<N<4 


N>5 




ADJ 








N> 1 


2-ASY 


N = 3,4 


N = 2,5 


6 < N < 14 


N> 15 



TABLE I: Catalogue of the four-loop zerology for 
SU(N), SO(N) and SP(2N) gauge theories with 
fermions transforming according to the fundamental 
and the 2-index representations. 



axis as an IR fixed point. The region where the new IR 
fixed point appears (on the negative coupling constant 
axis) ends before asymptotic freedom is lost. 



The defining feature shown in Fig. lc for topology C 
is that the appearance of two more merging points at 
negative values of a. 



In Fig. Id topology D, one observes that the IR zero at a 
positive value of the coupling reachers infinity at a finite 
value of the number of flavors, which is the distinctive 
feature of this topology. 

A new feature at the four-loop order is that two posi- 
tive nontrivial zeros, one IR and the other UV, can emerge 
simultaneously and can annihilate at a particular value 
of tif. At the two-loop level this feature does not exist 
and, in particular, no nontrivial ultraviolet fixed point is 
seen. 

As an example where these topologies arise we con- 
sider SU(N) with fundamental fermions as function of 
N. For N = 2 and 3 the topology A occurs. Increasing 
N the maximum value reached by the positive UV zero 
increases and for N = 4 it reaches infinity and there- 
fore it enters topology B. Increasing N further the local 
maximum of the IR negative zero-curve increases till it 
pinches the UV negative zero line for N = 11 entering 
topology C. Topology D is not realized in this case. On 
the other hand any SU(N) gauge theory with N > 2 
fermions and fermions in the adjoint representation lead 
to topology D. 

In table IT] we catalogue the four-loop zerology for 
SU(N), SO(N) and SP(2N) gauge theories with fermions 
transforming according to the fundamental and the 2- 
index representations. 



II. CONFORMAL WINDOW 

The conformal window is defined as the region in the- 
ory space, as function of number of flavors and colors 
where the underlying gauge theory displays large dis- 
tance conformality for a positive value of the coupling a. 
rif constitutes the upper boundary of the conformal win- 
dow and the lower boundary here is estimated by iden- 
tifying for which number of flavors the theory looses the 
infrared fixed point at a given number of colors. Due to 
the fact that we are using a truncated beta function the 
true window will be quantitatively different. 

We summarize the results for the SU(N) gauge groups 
in Fig.|2]for the fundamental, two-index symmetric, two- 
index antisymmetric and adjoint representation. The 
conformal window at the four-loop level is considerably 
wider, for any representation, when compared with the 
Schwinger-Dyson results 00 or the one obtained us- 
ing the critical number of flavors where the free energy 
changes sign, as suggested in JSJ. For completeness we 
also show the conformal window for the orthogonal and 
symplectic gauge groups respectively in Fig.|3]and Fig.l4] 
There is a universal trend towards the widening of the 
conformal regions with respect to earlier estimates using 
other nonperturbative methods. 



A. All-orders beta function comparison 

We have recently demonstrated the existence of a 
scheme in which the all-orders beta function is 0: 



jg(fl) _ a llC A -2T F w / (2 + A F y) 
3 \-2af x C A 



with 



A F = 1 + 



7_Ca 
11 C F ' 



(3) 



(4) 



The scheme independent analytical expression of the 
anomalous dimension of the mass at the IR positive zero 



is 1 : 



y 



llC A -AT F n f 



(5) 



We plot, for reference, in Figs. |2j |3j [4] the lines corre- 
sponding to this anomalous dimension equal to unity. 
These are the solid thick curves for the different repre- 
sentations. These lines could be viewed as the lower 
boundary of the conformal window if it is marked by 
the anomalous dimension to be unity. The size of these 
regions are consistent with the ones derived via gauge 



1 The same anomalous dimension was introduced in ] 10 1 motivated 
by the AdS/QCD correspondence. 



3 



CP 3 - Origins 



Particle Physics & Origin of Mass 




FIG. 2: Conformal window for SU(N) groups for the fundamental representation (upper light-blue), two-index 
antisymmetric (next to the highest light-green), two-index symmetric (third window from the top light-brown) 
and finally the adjoint representation (bottom light-pink). The lower boundary corresponds to the point where the 
infrared fixed point disappears at four loops. The solid thick lines correspond to the number of flavors for which 
the all-orders beta function predicts an anomalous dimension equal to unity. 
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FIG. 3: Conformal window for SO(N) groups for the fundamental representation (upper light-blue), two-index 
antisymmetric (which is the adjoint and second from the top (pink-region)), two-index symmetric (bottom window 

in light-brown). 



SP(2N) 
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FIG. 4: Conformal window for SP(2N) groups for the fundamental representation (upper light-blue), two-index 
antisymmetric (next to the highest light-green), two-index symmetric, i.e. the adjoint, (bottom window in light- 
pink). 
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FIG. 5: Anomalous dimension of the mass, at the in- 
frared fixed point, for SU(3) as function of the number 
of fundamental flavors at two loops (upper brown- 
curve), three-loops (second curve from the top in ma- 
genta), all-orders (dashed-curve in black), four-loops 
(bottom curve in blue). 



dualities in IITTl [T2l . Gauge duals have also been em- 
ployed in [14, 15J to compute important physical corre- 
lators such as the conformal S parameter |16|. 



B. Four-loop anomalous dimensions 

We plot, for illustration, the anomalous dimension of 
the mass for the SU(3) gauge theory, as function of the 
number of fundamental flavors, at the IR positive zero 
in Fig. [5] The three solid lines correspond respectively, 
from top to bottom, to the two-, three- and four-loop 
results. Perturbation theory is reliable only in a small 
range of flavors near rip A similar behavior is observed 
for any other gauge group, matter representation and 
different number of colors. We note that the perturbative 
analysis of the anomalous dimension appeared in IfTTl , 
while this paper was being completed. There it has also 
been noted that the anomalous dimension, to this order 
in perturbation theory, is smaller than for the three and 
two-loop case. 

Having at hand an all-order scheme-independent re- 
sult, we compare it with the perturbative one. The 
dashed line, in Fig. [5] is the all-order anomalous dimen- 
sion from |5j. It is striking that the all-order result is 
much more well behaved than the four-loop predictions 
which, in this example, reach large and negative values 
long before loosing the IR positive zero. 

Due to the phenomenological interest in models of 
minimal walking technicolor |6j [ISJ [191 we report the 
anomalous dimension at the fixed point also for the 
SU(2) gauge theory with two-adjoint fermions in Fig. [6] 
These theories are being subject to intensive numerical 
investigations via lattice simulations l20" W44l . 




0.4 0.5 0.6 0.7 0.8 0.9 1.0 



FIG. 6: Anomalous dimension of the mass, at the in- 
frared fixed point, for SU{2) as function of the num- 
ber of adjoint Dirac flavors at two loops (up green- 
curve), three-loops (second curve from the top), four- 
loops (third curve in blue), all-orders (dashed curve in 
black). 



III. ASYMPTOTIC SAFETY AT LARGE n f 

To the four-loop order a positive UV zero appears for 
a sufficiently large number of flavors. We have already 
observed that the value of the zero as function of number 
of flavors decreases monotonically as n~ 2 ^ 5 at four loops. 

In fact, it is possible to generalize this behavior to any fi- 
nite order in perturbation theory. Consider the equation 
for the zeros of the beta function in which the leading 
powers in the number of flavors are made explicit: 

CO 

b n f + Y^hn k f a k = , (6) 

k=l 

where b = fio/nf and hi = fik/n^. We used the fact 
that the first and second coefficient of the beta function 
are linear in the number of flavors and, in general, the 
successive coefficients have one extra power of rif ||45l . 
Therefore the coefficients b/ c are finite at large number of 
flavors. 
We define: 

x = rifa, (7) 
and the equation at any fixed perturbative order P reads: 

p 

b n f + Y J hx k = . (8) 

k=l 

At large n f the solution approaches: 
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There are P complex solutions for x lying on a circle in 
the complex plane. A positive solution exists only if bp is 
positive at large tif. This is indeed the case, at the four- 
loop order, for any gauge theory showing that the UV 
positive zero vanishes as n~ 2 ^ . If this UV zero persists to 
higher orders its location will change albeit will vanish 
faster as a function of tif when increasing P. 

Interestingly it is possible to sum exactly the perturba- 
tive infinite sum for the beta function, at large of number 
of flavors given that the leading coefficients are known. 
The result is: 



ArifT F a 2 



P(a) = 1 | H(x) | 



n f 



0(nf) . 



(10) 



The explicit form of H(x) can be found in |45|. The 
important feature, here, is that H(x) possess a negative 
singularity at x = 3n/Tp. This demonstrates that there 
always is a solution for the existence of a nontrivial UV 
fixed point at the leading order in rif for the following 
positive value of the coupling: 



«uv 



371 

T F n f 



(11) 



The function H(x) has also other singularities which 
might signal the presence of new zeros which we will 
not consider here, but that would be worth exploring. 
Higher order terms in n~ l can, in principle, modify 

the result if the singularity structure is such to remove 
or modify its location. 

A more complete discussion of the singularity struc- 
ture of the coefficients of the n~ l expansion has appeared 
in 1 45 1 also for QED. It seems plausible that the smallest 
UV fixed point is an all-orders feature. 



according to distinct representations of the underlying 
gauge group, at the four-loop order. The topology of the 
zeros of the perturbative beta function has been investi- 
gated. We discovered that only four distinct topologies 
are sufficient to classify the gauge dynamics of any the- 
ory. 

At the IR stable fixed point, for positive values of the 
a coupling, we computed the anomalous dimensions. 
We have also shown that these are well described by the 
all-orders beta function for any theory. 

Finally, by investigating the large n f limit we argued 
the possible existence, to all orders, of a nontrivial UV 
fixed point for any non-Abelian gauge theory at large 
number of flavors. 
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IV. CONCLUSIONS 

We presented the general features of the phase dia- 
gram for any gauge theory with fermions transforming 



Appendix A: Group factors and perturbative coefficients 



The four-loop beta function coefficients are 121: 



h 



11^ 4^ 

y C\ - 4C F T F n f - ^C A T F n f , 



2857 
54 



C\ + 2C 2 T F n f - 2 ^-C F C A T F n f 



(Al) 



1415 44^ ~ 2 158^ ^ , 

C 2 A T F n f+ -C F T 2 F n 2 f+ —C A T 2 F n 2 , 



W150653 44 \ riT I 



27 

39143 136 
3 



27 

Cb) + C 2 A C F T F n f {^ - ^Ca) + C A C 2 F T F n f (- 



81 



4204 352 
+ — 



, 7 /7930 224 \ 7 /1352 704 \ „ „ - W 17152 448 , \ 

+ 46C F T F n f + C 2 A T 2 F n}(— + — C 3 ) + C 2 F T 2 F n) (— - — C 3 ) + C A C F T 2 F n 2 f (— + — C 3 ) 

1232^-5 3 df d df d ( 80 704, \ df cd df d / 512 1664, \ ,df cd df cd ( 704 



243 



N A 



80 704 \ u F u7 /512 1664 \ 



N A 



27 
424 
243 ( 
512. 



+ —C A T 3 F n 3 f 



704 512 r \ 
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Here Cx is the Riemann zeta-function evaluated at x, T a F 
with a = 1, . . . ,Np are the generators for a generic repre- 
sentation F with dimension N F . The generators are nor- 
malized via tr(T F T b ) — Tp5 ab and the quadratic Casimirs 
are [T F T a F ]jj = C F bij. The subscript A refers to the ad- 
joint representation in the formulae in the text. Here the 
number of fermions is indicated by n f. 

The symbols d a F bcd are the forth-order group invariants 
expressed in terms of contractions between the following 
fully symmetrical tensors: 



jabcd 
F 



+T F T F T F T b F + T F T d F T b T F + T F T F T F T b ](A2) 



For readers' convenience we provide in table|Tl]the rele- 
vant group factors. 



The coefficients of the anomalous dimension to four- 
loops are J3l : 



Yo 

n 

72 
73 



3C F (A3) 

3 r2 97 10 
-C F + yC f C A - — C F T F n f 

™d - ™C\C A + ^C F Cl + CjT F n f (-46 + 48C 3 ) + C F C A T F n f (-^ - 48C 3 ) - ^Tjnj 

C t (-f- 336C 3 ) + C F C A (If? + 316C 3 ) + C F Cl (- - 152C 3 + 440C 5 ) + + ^C 3 - 440t 5 ) 

280 „ \ m / 8819 



(/8f) \ / 8819 \ 

-— + 552C 3 - 480C 5 ) + C f C a T fM/ — + 368£ 3 - 264£ 4 + 80C 5 j 

+C F C 2 A T F n f (-^--^3 + 264C 4 + 400C 5 ) + CjTjnj (f - 160C 3 + 96C 4 ) 

^- + 160C 3 - 96C 4 ) + C F T%n) \-~ + -^C 3 ) + (-32 + 240C 3 ) + n f F — F (64 - 480Q 



jabcd rfabcd 



r 



The results of \A1\ and I A3 1 ares valid for an arbitrary 
semi-simple compact Lie group. The result for QED (i.e. 
the group U(l)) is included in Eq. ||A3} by substituting 



C A = 0, d abcA = 0, C F = 1, T F = 1, (df ci f = 1, N F = 1. 

These coefficients were obtained in an arbitrary covari- 
ant gauge for the gluon field and are gauge independent. 
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Rep. 



Nr 



df cd df d /N F 



df cd df cd /N F 



SU(N) 



(N-1)(N+1)(N 4 -6N 2 +18) 
96N 3 

^(1^ + 36) 

(IV-l)(N+2)(N 5 +14N 4 +72N 3 -48N 2 -288N+576) 
48N 3 

(N-2)(N+l)(N 5 -14N 4 +72N 3 +48N 2 -288N-576) 
48N 3 



FUND 
AD} 

2-SYM 
2-ASY 



N 
N 2 -l 

\N(N + 1) 
±(N-1)N 



1 




N 2 -l 




2 




2N 




N 




N 




ZV+2 
2 


N 


-i + 


1 


N-2 
2 


N 


2 
N 


1 



i (N _l )(N+ l)( N 2 + 6) 

^N 2 (N 2 + 36) 
i(N - 1)(N + 2) (N 2 + 6N + 24) 
g(N - 2)(N + 1) (N 2 - 6N + 24) 

SO(N) 



FUND 
AD} 
2-SYM 


N 

±(N-1)N 
\ (N 2 + N - 2) 


1 

N-2 
N + 2 


N-1 
2 

N-2 
N 


(N-1)(N 5 -10N 4 +315N 3 -1250N 2 +1840N-1408) 


(N-1)(N 4 -2N 3 +107N 2 -106N+128) 


384(N 2 -N+4) 
N 6 -18N 5 +875N 4 -6170N 3 +17600N 2 -25728ZV+18944 


384(N 2 -N+4) 
N 6 -18N 5 +875N 4 -6170N 3 +17600N 2 -25728N+18944 


192(N 2 -N+4) 
N(n 5 -4N 4 +723N 3 -2576N 2 +2752N-1792) 


192(N 2 -N+4) 
N(jV 5 +12N 4 +787N 3 +1824N 2 +1344JV+3328) 


192(N 2 -N+4) 


192(N 2 -N+4) 


SP(2N) 


FUND 
AD} 
2-ASY 


2N 
N(2N + 1) 
N(2N - 1) - 1 


i 

2 

N + l 
N-1 


|(2N + 1) 
N + l 
N 


(2N+l)(4N 5 +20N 4 +15N 3 -50N 2 +10N+10l) 


(2N+l)(8N 4 +8N 3 -26N 2 -14N+49) 


96(2N 2 +N+2) 
4N 6 +36IV 5 +125N 4 +85N 3 -100N 2 +216N+434 


384(2IV 2 +N+2) 
4N 6 +36N 5 +125N 4 +85IV 3 -100N 2 +216N+434 


24(2N 2 +N+2) 
N(4N 5 +8IV 4 -27IV 3 -62N 2 +238N+374) 


24(2N 2 +N+2) 
N(4N 5 -24N 4 +37N 3 +138N 2 +6N-386) 


24(2N 2 +N+2) 


24(2N 2 +N+2) 



TABLE II: Relevant group factors for SU(N), SO(N) and SP(2N). 
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